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INTRODUCTION 



I . 

Performance prediction codes for solid rocket motors 
model two phase flow losses as functions of particle size. 

In addition, particle size within the grain port strongly 
affects the damping of combustion pressure oscillations. At 
present these models are based on particle size data from 
collected exhaust samples [Ref. 1] . However, particle size 
varies with position in the motor and otner parameters 
(pressure, propellant formulation, nozzle design, etc.) . 
Therefore, experiments to determine how particle size varies 
in the actual flow environment of the motor (i.e., across the 
nozzle) are needed to validate the models for two phase flow 
losses. Cramer [Ref. 2] and Karagounis [Ref. 3] provide a 
good summary of the subject and the Naval Postgraduate School 
Combustion Laboratory effort to obtain particle size data 
across the exhaust nozzle of a solid rocket motor. 

The method used in this continuing effort was the dif- 
fractively scattered light technique. The diffraction patterns 
of light scattered by particles are analyzed to determine the 
volume to surface mean diameter [Refs. 4 through 11] . This 
method has the disadvantages that size distributions cannot 
be easily determined and particles larger than some threshold 
size will not be detected due to the exceedingly small angles 
as which they scatter light. However, it has the advantage 
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that it is r.on-intrusivs anc m z-.eory can be used m tue 
internal motor environment. 

Use of this method for particle sizing in solid rocket 
motors at the Naval Postgraduate School Combustion Laboratory 
was begun by Karagounis [Ref. 3] . fhe apparatus was sub- 
sequently redesigned and the data acquisition equipment up- 
graded with the introduction of the Hewlett Packard 3054A 
data acquisition system with an HP 85 as the controlling 
computer [Ref. 12] . The investigation by Cramer and Hansen 
followed and showed that propellent composition can limit 
the application of the technique. Large particulate combustion 
products in the flow made particle size data difficult to 
obtain. This was especially true if only one measurement of 
the scattering profile was made during a test firing. 

To address this problem in the present study several 
improvements were made. A cleaner burning propellant was 
obtained to reduce char agglomerates in the exhaust products. 

A more statistically valid data sample (multiple measurements 
during a single test) was made possible with added momory 
in the data acquisition equipment and a pacing circuit which 
allowed full use of this memory. Data reduction was also 
improved with a Hewlett Packard 9836S computer combined with 
a more recently developed approach to particle sizing 
presented by Buchele [Ref. 13] . This method is discussed 
later in depth. 
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The focus of tnis cr.es is crcjacc was the following: 

(1) Implement the HP 9836S as the system controller. 

(2) Expand the multiprogrammer memory in order to obtain up 
to eight consecutive scans of the diode arrays during 

a test firing. 

(3) Improve data reduction techniques by the method of 
Buchele [ Ref . 13 ] . 

(4) Modify the apparatus and experimental procedures to 
improve the angular resolution and to reduce extraneous 
light . 

(5) Certify the proper functioning and accuracy of the 
apparatus prior to actual motor testing. 
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THEORETICAL BACKGROUND 



II . 

A. GENERAL DISCUSSION 

The completely general theory of scattering was developed 
by Mie and is presented by Van de Hulst [Ref. 14] . The 
light scattering characteristics for spherical particles of 
any size are fully described in a mathematical format. The 
Mie scattering functions contain Legendre polynomials and 
spherical Bessel functions and fully treat the phenomena of 
reflection, refraction, diffraction, and extinction. The 
full theory is most often applied when particle size is 
approximately the same as the wavelength of the incident light. 
Van de Hulst [Ref. 14] calls this the regime of Anomalous 
Scattering . 

For particle sizes much smaller than the wavelength of 
light the Mie equations simplify to a form which is more 
dependent on the index of refraction of the particles and 
less dependent on particle size. This is called Rayleigh 
Scattering . 

The study of particle size behavior in solid propellant 
rockets mainly covers sizes much greater than the wavelength 
of light. Scattering by large particles such as these is 
described adequately by Fraunhofer diffraction. 
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B. APPLICATION TO LARGER PARTICLES 

The ringed diffraction pattern generated by a hole in 
a mask, or a number of particles of the same size is de- 
scribed by the equation: 

2J ( ot 6 ) 2 

1(6) = [ ] 

a0 



where : 

1(6) describes the relative intensity of the scattered 

light at an angle theta (0) 

J^(a0) is the Bessel function of the first kind 
TT D 

a = t — is the particle size parameter for diameter 

D and wavelength of light lambda (X) . 

Measuring the particle size for a monodispersion can be 
accomplished by measuring the angular position of a dark or 
bright ring in the diffraction pattern. For a dark ring the 
zero of the Bessel function corresponding to the ring is set 
equal to (a0) and particle size is determined directly [Ref. 5] 
For bright rings one sets (a0) equal to the corresponding 
maximum of the Bessel function and solves for the diameter. 

The above method is not used for polydispersions since 
the discrete rings are not observed. However, Dobbins, et al . 
[Ref. 5] introduced a significant improvement in the dif- 
fractively scattered light method of particle sizing. They 
found that although the method was not directly able to 
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determine distributions of sizes, the volume to surface 
mean diameter defined by 



(D)D 3 dD 
0 r 

/^°°N (D)D 2 dD 
0 r 



where : 

N (D) is a distribution function describing the pro- 
portion of particles with diameter (D) in the sample, 
could be accurately measured. 

A curve for sizing polydispersions was presented which was 
used by Cramer and Hansen [Refs. 2, 12] . 

Two phase flow losses are often calculated in terms of 
If the distribution of sizes in the polydispersion is 
well behaved then and other diameters such as can 

be easily related [Ref. 15] . Reference 5 reported that very 
small particles in the distribution have a minor influence 
on the scattering profile. This makes the measurement 
technique promising for the two phase flow loss study 
since very small particles do not contribute significantly 
to these losses and so are of less interest. 

Roberts and Webb [Ref. 6] essentially confirmed the 
conclusions of Reference 5 and presented a similar curve for 
use in sizing. 

More recently, Buchele [Ref. 13] gives a good summary of 
experimental techniques for particle sizing by measuring 
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diffracted light. One point of interest: in ms report is 
that he represents the scattering profile of a polydispersion 
with a function which closely approximates the curves of 
References 5 and 6. 

I (6) = EXP - ( . 57a6 ) 2 

n 

This function from Reference 13 and the curve from 
Reference 5 were both used in the present study to evaluate 
the apparatus to be used with solid propellant rocket motors. 

An additional detail of measuring scattering profiles is 
covered by Van de Hulst [Ref. 14] . The wavelength of light 
used in the scattering calculations depends on the index of 
refraction of the medium containing the particles. The 
wavelength used in all calculations must be: 



where : 

Aq is the wavelength of light in a vacuum and, 

M is the refractive index of the medium with respect 
to a vacuum. 

Thus, the size parameter (a) becomes: 




and the beam spread parameter becomes: 
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Another consideration is presented by Gumprecht and 
Sliepcevich [Ref. 4] . Light scattered by particles in a 
medium is refracted as it crosses each interface of the 
container holding the medium. This is discussed further in 
the section on calibration and evaluation of the apparatus . 

Additional complications arise with the full treatment 
of the index of refraction of the particles with respect to 
the medium. But, for Fraunhofer diffraction alone this aspect 
can be neglected. 

C. RESTRICTIONS AND SOURCES OF ERROR 

Some restrictions on use of the method are described by 
Dobbins, et al . [Ref. 5] and were satisfied as described by 
Cramer [Ref. 2] . These are related to the size of particles, 
the distance to the detector, and some phenomenon covered in 
the rigorous Mie theory. 

One must keep in mind also that the curves developed for 
polydispersions are based on the Upper Limit Distribution 
Function of Mugele and Evans [Ref. 7] . This means that no 
particles with size greater than approximately ten times the 
mean should be in the sample [Ref. 13] . This appears to be 
a mild restriction. Van de Hulst [Ref. 14] describes the 
criteria for single scattering and a simple test to verify 
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it. In general, as long as the scattered intensity is 
proportional to the number of particles the mathematics 
remain simple. 

Sources of error of the dif f ractivelv scattered light 
method are covered by Buchele [Ref. 13] and are presented 
here . 

(1) Inaccuracy of angular measurement or the limited ability 
to resolve small angles and, 

(2) Inaccuracy of the intensity measurement due to extraneous 
light . 

Extraneous light includes all light other than scattered 
light from the particles-. Some examples are scattering from 
an aperture or dirty test section windows. Refraction of the 
beam due to gas density gradients and image point broadening 
from turbulence are others. Laser speckle is also extraneous 
light . 

The sources of error addressed in this investigation are 
discussed in the- related portions of the paper. 
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Ill . 



EXPERIMENTAL APPARATUS 



A photograph of the apparatus is presented in Figure 1. 

A schematic is presented in Figure 2. The light scattering 
equipment was mounted on two optical benches. Components for 
measurements in the exhaust plane were mounted on one bench. 
The other bench held the equipment associated with the motor 
cavity. The light source was an eight (8) milliwatt Helium 
Neon laser mounted on the exhaust bench. A collimated beam 
was required so a spatial filter/collimator was used. A 
modification to. this collimator is discussed later in this 
section. The collimated beam passed through a cube beam 
splitter and the second beam was diverted to a 90 degree 
prism on the other bench. The original beam continued through 
the motor exhaust plane. The other beam was routed through 
the nitrogen-purged glass windows in the motor housing. 

Each beam was then intercepted by a physical stop located 
in front of its set of receiving optics. The further the 
stop was placed from the test section, the smaller the angle 
at which scattered light could be measured. In this 
apparatus, the stops were placed approximately 30.5 
centimeters from the exit plane of the test section. This 
allowed a minimum angle of approximately .008 radians to be 
measured. Light scattered at angles greater than this was 
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The stop served to keep the transmitted beam out of the 
measuring optics and thus reduce extraneous light. The 
stops also improved optics alignment. This is discussed 
under calibration and evaluation of the apparatus. 

The scattered light passed through a narrow pass filter 
which admitted only light of the Helium Neon frequency. 

This filter served to reduce extraneous light from the external 
surroundings . 

An objective lens of 50 centimeter focal length was 
located behind the narrow pass filter. This lens imaged 
onto a photodiode array the scattering profile of the particles 
in the test section. The shadow of the beam stop was also 
imaged since the stop was between the test section and the 
objective lens. This was a limitation which is discussed 
under calibration and evaluation. 

The photodiode arrays were the same units used by Cramer 
and 'Hansen [Refs. 2,12]. Each array contained 1024 silicon 
photodiodes on a single chip with 25 micron spacing. The 
accompanying circuits provided a sampled and held output 
which was essentially analog except for switching transients. 

At the end of each diode scan there was a delay before the next 
scan. During this delay the diodes were reset and allowed 
to measure the intensity of the scattering profile again. 

The scanning of the diode array repeated continuously. The 



actual sampling time of the array was about 34 milliseconds 
with a delay between scans of about 6 milliseconds. 

The 50 centimeter focul length of the objective lens 
combined with the dimensions of the diode array provided a 
half angle field of view of about 3 degrees for mediums of 
refractive index near one. The effective field of view was 
reduced to about 2.3 degrees for calibrations when the 
refractive index of a Plexiglas container and water was 
taken into account. 

The laser beam collimator mentioned previously at 
first produced a beam one centimeter in diameter. A lens 
in the collimator was changed to reduce the beam diameter for 
several reasons. Extraneous light would be generated if 
a large beam impinged on the aperture of the motor test 
section window. Also, if the aft beam was larger than the 
motor exhaust jet it would be refracted in the density 
(and refractive index) gradient between the exhaust and air. 

The last part of the apparatus was the rocket motor 
itself. It was the one used by Cramer and Hansen [Refs. 2,12] 
and in the present study served only for aligning the optics. 
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IV. 



DATA ACQUISITION SYSTEM 



A. NEW CONTROLLER 

Hansen [Ref. 12] describes the major components of the 
Hewlett Packard 3054A data acquisition system. A list of the 
manuals relevant to this study is in Table 1. The HP 85 
computer used by Hansen and Cramer was replaced with an 
HP 9836S as system controller. This newer computer has far 
more capability than the HP 85, including a choice of more 
powerful operating systems. The system used for this 
study was Basic Extended 2.1. 

The data acquisition program written by Hansen needed 
minor modification to acquire multiple consecutive scans of 
the photodiodes. Some different I/O commands such as those 
which transfer data to the disk were also incorporated. 

The revised version of this program is listed in Appendix A. 

A general flow chart is presented in Figure 3. 

The 9836S has two internal disk drives which were used 
to store the data after acquisition. The data from both 
diode arrays was stored in the same file. The eight (8) 
scans of the motor cavity were first, followed by the four (4) 
exhaust scans. 
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B. MODIFICATIONS 

The memory capacity of the Multiprogrammer unit was 
increased so that multiple consecutive scans of the diode 
arrays could be recorded during a motor firing. This would 
provide a more statistically valid measurement of particle 
size. Fluctuations of scattered light intensity for a 
polydispersion need to be integrated over time or averaged 
to provide a more appropriate measurement. 

In order to fully use the memory added and make data 
management easier the data acquisition system needed to be 
modified. The memory consisted of three (3) cards, each 
with a capacity of 4096 values. The fact that this was a 
multiple of 1024 (the exact number of photodiodes) meant 
that the idle period between scans needed to be excluded 
from the data. If this was not done, one (1) less scan per 
card would have been acquired and locating the scans in 
the overall block of memory would have been more difficult. 

It was also necessary to chain two of the memory cards 
together in a way which would allow one card to be filled 
and then the other. A schematic of the data acquisition 
system is presented in Figure 4. 

C. PACING AND MEMORY CONTROL CIRCUITS 

The timing clock and blanking pulse of the photodiode 
circuitry provi-ded the means for pacing data acquisition. 
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Specific results of the modif icaticr.s were: 

(1) Memory space was fully utilized and management of the 
multiple scans made easier. 

(2) A/D conversions of the data were made exactly when a 
diode's output was on line and steady. Thus, the analog 
filter used in the previous study was no longer 

needed to suppress the switching spikes on the data line. 

The following is a description of the signals and circuits 
used to modify the data acquisition system. All voltage 
levels were TTL . A timing diagram in Figure 5 shows the 
relations between signals. A schematic of the circuit is 
presented in Figure 6 . 

The clock pulse was a positive going spike at a frequency 
of about 30 KHz. This clock controlled all circuits of 
the photodiodes. It ran continuously, even during the blank 
period between scans when the diode output was clamped at 
zero volts. 

The blanking pulse was a signal which fell to zero at 
the beginning of each scan. It then went positive at the 
end of the scan and remained high until the next scan began. 

The clock pulse was used to drive a pulse shaper (mono- 
stable multivibrator) . This ensured that the voltage levels 
through the rest of the circuit would not accidentally fall 
below the TTL threshold. The pulse width of the shaper 
was adjusted so . that the negative going edge of each pulse 
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would occur after the switching transient on the data line 
had decayed. This negative going edge would eventually 
trigger the A/D converter to store the output of each diode. 

The blanking pulse was inverted and connected to an AND 
gate along with the pulse shaper output. The output of this 
AND gate is shown in Figure 5 as the pulse shaper signal held 
low between scans of the diode array. This was the basic 
signal which paced data acquisition. 

This basic trigger signal was connected to an AND gate 
along with the output of the Multiprogrammers Timer Pacer 
card. In this way the trigger would not reach the A/D until 
the Timer Pacer output a pulse. This enabling pulse from 
the Timer Pacer was at least as long as the time for eight 
scans of a diode array. 

The controller programed the Timer Pacer to produce the 
pulse when the Timer Pacer received a trigger from the blanking 
pulse. In this way, data acquisition began at the start of 
a scan and no data was taken during the time between scans . 

The circuit to chain the memory cards together was 
basically an OR gate used as a negative logic AND gate. The 
end of conversion (EOC) signal of the A/D and the (FULL) signal 
of one memory card were connected to the gate. When both 
signals went low the second card was then able to store data 
from the A/D. This arrangement is shown in the schematic of 
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the data acquisition system in Figure 4. The automatic lock- 
out feature of the memory cards when full, and the relatively 
slow rate of data throughput made it unnecessary to control 
other handshake lines [memory card manual. Table 1] . 

The circuit was designed to handle four (4) memory cards, 



so no 


modification will be 


necessary if one 


more card i 


added 


to 


the 


system . 


This 


would provide an 


additional 


scans 


of 


the 


exhaust 


beam . 
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V. DATA REDUCTION 



The data reduction programs written by Hansen were not 
used for this study. The new computer lended itself to 
another approach. Its memory capacity made it unnecessary to 
chain programs together and polynomial curve fitting was 
eliminated in favor of interactive graphics. Avoiding 
polynomial fits preserved the nature of the raw data so that 
one had a better feel for the parameters. The data reduction 
program "RDC" is listed in Appendix A. Figure 7 is a general 
flow chart for the program. 

The following is a description of the program. The user 
was first prompted for values needed to analyze a given 
data set. For example, the wavelength of the laser used and 
the index of refraction of the medium must be known for 
any data set. Next, one had the choice of reducing raw 
data scans or reviewing a reduced data file. For raw 
data one chose either the exhaust or motor cavity 
beam data. 

Raw data was plotted on the CRT and any obviously erroneous 
scans were excluded from further reduction. The valid scans 
were averaged to obtain a mean scattering profile. The mean 
intensity profile taken before particles were introduced 
was then subtracted from that taken, with particles present. 
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This corrected for the characteristics of individual photo- 
diodes and extraneous light which was independent of the 
particles . 

A symmetric moving-average-type of digital filter was 
then applied to the profile to achieve some smoothing. 

This type of digital filter was chosen for simplicity and 
because it does not have the phase lag of analog filters 
[Ref. 16] . Preserving the phase of the data was necessary 
to retain angular resolution. Another advantage of filtering 
in the software rather than hardware was that raw data 
files remained unmolested. 

The scattering profile was then analyzed using interactive 
graphics. If earlier, one chosen to review a reduced file, 
program execution began here. 

One had to normalize a scattering profile in order to 
compare it to the theoretical curves for polydispersions. 

The scattered intensity on the centerline of the beam was 

the correct value to use for normalization but was unmeasurable 

due to the beam's presence. 

The other unknown was, of course, the particle size. These 
two variables (centerline intensity for the measured profile 
and D 32 f° r the theoretical profile of normalized intensity 
vs. ( 6 ^)) were adjusted using interactive graphics until 
the curve for polydispersions coincided with the data. In 
this way the mean diameter of particles was determined. 
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The second reduction technique used was the direct 
application of the method presented by Buchele [Ref. 13] . 

The equation for the polydispersion curve: 

I (9) = EXP - ( . 57 a0 ) 2 
n 

was applied at two points of the scattering profile. This 
gave : 

12/11 = EXP -D 2 [ ( © 2 " 0^ )( ' 57 VA) 2 ] 

Solving this for the diameter gave: 

D = [-L (12/11) (A/.57") 2 /(6 2 - 0 2 )p 
n 2 1 

The computer would sweep through the data using many 
values of thetal (0^) along with several angle ratios to 
determine theta2 (© 2 ) • The results were presented graphically 
as particle size vs. thetal for each angle ratio (Q^/B ) . 

In actual practice the range of useable angles depends on 
the appratus, and the quality of the data. Therefore, in order 
to interpret the results one must have previously inspected 
the data. The interactive graphics routine was well suited 
to this and provided a hard copy for inspection. 

After reducing a set of raw data the mean scattering 
profile was stored on disk for later review. 
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VI . 



CALIBRATION AND EVALUATION 



A. IMPROVEMENTS 

The geometry of the apparatus used in the investigation 
by Cramer and Hansen is compared with that of this study 
in Figure 8. In the previous study the transmitted beam 
was allowed to enter the receiving optics. The beam was 
focused off the diode array a few millimeters from the first 
diode. This was necessary to avoid damaging the diodes 
but introduced some uncertainty in angle measurements. 

The intense image of the beam along with scattered light from 
the receiving optics produced a high level of extraneous 
light. In the present study, stops were used to intercept 
the beam before reaching the receiving optics. These stops 
provided several advantages. A high intensity beam could 
be used while producing little extraneous light. Optics 
alignment was also improved. This reduced error in angle 
measurement. Alignment was accomplished using a neutral 
density filter to reduce beam intensity and protect the 
diodes. A schematic of the apparatus is in Figure 2. The 
laser, collimator, beam splitter and prism were positioned so 
that the beams passed through the appropriate measurement 
areas. The narrow pass filter and imaging lens were then 
positioned so that the beam entered on the centerline. The 
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photodiode array was then moved using a three-axis micrometer 
so that the focused beam fell on the first diode. The beam 
stops were then put in place and the neutral density filter 
removed. In this way measurements commenced exactly from 
the optical axis of the beam. 

Procedures were also refined to account for the bending 
of light rays as they passed through the walls of the particle 
container. As noted earlier, the index of refraction of the 
container and the medium containing particles affects 
scattering measurements. A Plexiglas box held the particle 
samples and a magnetic stirrer kept the samples suspended 
in water. The index of refraction of the Plexiglas and 
water combination was measured using a simple technique. 

A microscope was used to measure the ratio of actual depth 
to apparent depth for Plexiglas and water. The index of 
refraction was determined to be 1.39. This value was 
applied to the data to convert the measured scattering profile 
to that actually produced in the medium containing the 
particles . 

B. RESULTS 

Calibration results are summarized in Table 2. Initial 
tests were done with two samples wihch were basically mono- 
dispersions of large particles. Figure 9 shows the measured 
profile of scattered light for glass spheres ranging from 
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37 to 44 microns in diameter. This profile was obtained by 
placing the focussed beam just far enough from the first 
diode to avoid saturation with no particles present. The 
diodes located at angles less than about .01 radians 
saturated. The first bright ring for particles of about 40 
micron diameter was visible near .02 radians. Figure 10 
shows a profile for the same particles, illustrating use of 
the beam stops to avoid diode saturation and improve angle 
measurements. In this case, the first diode was located 
exactly on the centerline of the beam as discussed above. 

Results for a sample of 53 to 63 micron glass spheres are 
shown in Figure 11. The center lobe was nearly completely 
missed but the first two bright rings were seen near .014 
and .022 radians. The first two dark rings near .01 and 
.019 radians were also seen. The method described earlier 
of setting the beam spread parameter ( ot 0 ) equal to the 
zeroes and maximums of J^(a0) was used to calculate a size 
of about 58 microns. Also shown in Figure 11 is the 
theoretical profile for a polvdispersion with D^-, = 54 
microns . 

Various polydispersions of either glass spheres or 
aluminum oxide powder were then tested. These polydispersions 
consisted of fairly large particles. Results are shown in 
Figures 12 through 17. These tests showed that the apparatus 
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had two distinct modes of operation. If the particle concen- 
tration was very high, or if large particles dominated the 
polydispersion, many of the diodes at the smaller angles 
would saturate. This left only the data at larger angles 
useable. When many diodes saturated, the theoretical curve 
given by Dobbins, et al . , was used to determine size. This 
was done because this curve was valid for the larger angles 
and lower relative intensities. The curve from Buchele 
[Ref. 13] was not valid for values of the beam spread 
parameter greater than three (3) . 

For low particle concentrations and/or small particles the 
data proved more accurate at the smaller angles. If no 
diodes were seen to saturate then one knew the measurement 
was in the higher intensity part of the center lobe. Here 
the curve given by Buchele was quite satisfactory for sizing. 

The smallest particles tested were five, ten, and twenty 
micron polystyrene spheres. The bright rings for these 
particles occurred at angles too large for the apparatus to 
measure. For these samples the diodes did not saturate. 

Both the Gaussian curve fit and the two angle method were 
used to obtain These results were especially consistent. 

It should be noted that the two-angle method uses the equation 
for the Gaussian. If the measured profile matched the 
Gaussian exactly, then would be the same for any (0^) 

and angle -ratio ( 82 / 02 ) employed. Some variations in 
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calculated <tue to the imperfect fit are obvious in 

Figures 19, 21 and 23. 

A scanning electron microscope was used to photograph 
the types of particles tested. These photographs are shown 
in Figures 24 through 27. Equation (1) was used along with 
these photographs to calculate some of the values of 
in Table II. Calculations of for the polystyrene were 

arrived at using the manufacturers data on size distributions. 
The photographs generally confirm the validity of the 
technique . 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



The results of the calibration tests showed that the 
apparatus is capable of accurately measuring mean particle 
size for a broad range of mean diameters. It was found that 
the technique was most accurate if the theoretical profile 
fit or the two-angle method were applied at the smallest 
possible scattering angles. 

The rocket exhaust is likely to attenuate the beam 
somewhat, reducing the problems related to diode saturation 
at small scattering angles. Thus, measurements should be 
possible using the high intensity part of the center lobe. 
This should make data reduction less ambiguous. Actual 
testing should begin with measurements at the exhaust 
plane of the motor. These should be compared with collected 
exhaust samples to validate the use of the apparatus in 
an actual motor environment. Measurements in the motor 
cavity would then be interpreted based on the correlation 
between exhaust samples and exhaust measurements. 

It is also recommended that the index of refraction of 
the combustion gases be investigated. A literature search 
for an estimate of the index of refraction would probably 
be satisfactory. 
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TABLE I 



ELECTRONICS MANUALS 

1. HP Memory Cards Model 6970B Operating Manual 

2. HP Timer Pacer Card Model 69737A Operating Manual 

3. HP Analog to Digital Converter Card Model 69736A 
Operating Manual 

4. HP Users Guide, "Using the 9826 and 9836 Computers with 
the 6942A Multiprogrammer" 

5. Basic Language Reference Guide with Extensions 2.0 for 
Series 200 Computers 
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TABLE II. CALIBRATION RESULTS 



Error 


microns 




LO 


rO 

CX] 


VO 


CXI 


p 

o 

p 

o 


3 to 5 


1 

i 

l 

1 

l 

1 


l 

I 

l 

l 

l 

l 


P 






















so e 






















£ CD 














CO 


o 






* H £ 




CO 










LO 


rO 






P (D 




C 


















CD P 


CXI 


O 


LO 


cn 






o 


O 






P 3 


rO 


P 


• 


• 






p 


P 






p (/) Q 


u 


p 




!p 


o 






CO 


LO 


cd cd 




H 






CXI 


p 


p 


CO 


CXI 


P 


U CD 




£ 










LO 


CXI 






oo 2 






















s 






















rP 






















w CD 






•K 


•K 


*K 












P 




CO 


* 


4C 




•K 


•K 


* 


P 


XT 


G cd 




C 














CX 


CX 


O t— 1 


CX] 


o 




cxi 


xD 








<D 


0) 


*H £ 


to 


p 


. 


. 


• 


. 


. 


. 


0) • 


CD • 


+J O Q 


u 


p 


o 


tH 


CO 


P 


LO 


CO CO 


CO CO 


cd * — i 




• p 




!— 1 


CXI 


r O 


LO 


CXI 


•H 


P 


G cd 




£ 














PP 


fP 


cru 






















w 






















CD 








xD 


o 


P 


rO 


r- 






t— 1 


CO 




xD 


!— 1 


rO 


p 


xD 


rO 






U CD 


G 
















uo 


O 


*H N 


o 




o 


o 


O 


O 


O 


o 


CXI 


uo 


P -H 


p 




p 


p 


P 


P 


P 


p 






P 00 


u 
















*1 


l\ 


cd 


• H 




rO 


xD 


LO 




rO 


!p 






CP 


£ 








iP 


rO 


LO 














CD 


CD 


CD 












CD 


t— 1 1 




G 


c 


G 








£ 


£ 


r— 1 


cd 




CD 


CD 


CD 








G 


G 


u 


• H 




P 


P 


P 


CO 


CO 


CO 


G CD 


G CD 


• rH 


P 




X 


X 


X 


CO 


CO 


CO 


•H ^ 


•H T3 


P 


CD 




P 


p 


P 


cd 


cd 


cd 


£ *H 


£ *H 


P 


P 




CO 


CO 


CO 


?p 


rP 


fP 


G X 


G X 


cd 


cd 




X 


X 


X 


O 


00 


O 


i — i O 


rP O 


Co 


S 




tP 


»P 


»P 








< 


< 








o 


o 


o 


















CP 


CP 


CP 













36 



Prom SEM Photos ** Prom Manufacturers Data 





Figure 1. 



Photographs of 



Light Scattering Apparatus 



37 






Hel ium - 

Neutral Neon 




U 

o 



33 



Figure 2. Schematic of Light Scattering Apparatus. 
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Figure 3. Flow Chart for Program "ACQDTA" . 
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Figure 4. Schematic of Data Acquisition System. 
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Figure 5. Timing Diagram for Data Acquisition. 




Figure 7. (Continued) Flow Chart for Program "RDC " . 
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Figure 7. Flow Chart for Program " RDC " . 
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Figure 7 . (Continued) Flow Chart for Program " RDC " . 
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Figure 8. Comparison of Two Geometries for Light Scattering Measurements. 
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Figure 9. 37-44 Micron Glass Spheres Using No Beam Stop. 
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Figure 10. 37-44 Micron Glass Spheres Using Beam Stop. 
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Figure 11. 53-63 Micron Glass Spheres Showing Diffraction Rings. 



CURVE FIT RESULTS 
INTENSITY vs.THETR 




C3GDCSCDGD£SGDC3CSCD£SC3 



CD I 

AXISN31NI aszinuwuoN 



0 

•H 

4 J 

03 

u 

4 J 

C 

a 

o 

c 

o 

o 

x: 

C 7i 
•H 



CN 



( 1 ) 

U 

P 

tn 

•H 



49 



THETR (rad) 



CURVE FIT RESULTS 
INTENSITY vs.THETR 




A1ISN31NI QSZinUNUON 



50 



Figure 13. 1-37 Micron Glass Spheres, Low Concentration. 
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50 Micron Aluminum Oxide, High Concentration. 
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Figure 15. 50 Micron Aluminum Oxide, Low Concentration. 
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Figure 16. 25 Micron Aluminum Oxide, High Concentration. 
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Figure 17. 25 Micron Aluminum Oxide, Low Concentration. 
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Figure 18. 5 Micron Polystyrene, Curve Fit. 
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Figure 19. 5 Micron Polystyrene, Two Angle Method [Ref. 13] . 
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Figure 20. 10 Micron Polystyrene, Curve Fit. 
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Figure 21. 10 Micron Polystyrene, Two Angle Method iRef. 13]. 
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Figure 22. 20 Micron Polystyrene, Curve Fit. 
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Figure 23. 20 Micron Polystyrene, Two Angle Method [Ref. 13] . 
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Figure 26. SEM Photographs of Glass and Polystyrene Spheres. 
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Figure 27. SEM Photographs of Polystyrene Spheres. 




APPENDIX A 



PROGRAM LISTINGS 



10 !********** AC9DTA ; FOR 9336S ****************** 

20 !********** ACQUIRES MULTIPLE SCANS GF BOTH DIODE ARRAYS ***************** 



30 !********** AND STORES THEN ON DISK ***************** 

40 !********** BERT HANSEN, KEN GRAHAM, KELLY HARRIS 1984 ***************** 

50 OPTION BASE 1 

SO CON QH112] ,Q2$112] ,T1$120] ,D1$112] ,D2$l 12] ,Address(4) ,E(4096) 

70 COM AG024) BUFFER ,B(409S) BUFFER, C(40?6) BUFFER ,D(4096) BUFFER 
80 OUTPUT 709; ’AR* ! ANALOG RESET 

90 ASSIGN gHultilO TO 72310 (CLEARS THE UAXE-UP SERVICE REQUEST 

100 ENTER SHulUlBjQql ,Qq2,0q3,Qq4 ) Qq5 ) Qq6 !0F THE MULTIPROGRAMMER 



110 


HAT Address= 


(0) 




120 


ENABLE INTR 7 




130 


CLEAR 722 






140 


PRINT USING 


■3/" 




150 

)* 

160 


PRINT "ENTER THE FILENAMES OF THE DATA FILES TO BE CREATED (e.q, RAU1 ,RAW2 


PRINT " 




12 CHARACTERS MAXIMUM, EACH 1 


170 


PRINT USING 






180 


PRINT * 


AN EMPTY 


DISK MUST BE IN THE LEFT DISK DRIVE 1 


190 


PRINT '■ 






200 


PRINT * 


THE DATA FILES WILL NEARLY FILL A DISK" 


210 


PRINT USING 


7" 




220 


HASS STORAGE IS " : INTERNAL , 4 , 


,0" ! CHANGE THIS LATER IF NECESSARY 


230 


INPUT * 


INPUT FILE NAMES NOW - (FILENAME! .FILENAMES) M ,D1$,D2$ 


240 


H5=1024 






250 


Zz$= INTERNAL, 4, 1" 


! STRING INDICATES MASS STORAGE 


260 


CREATE BDAT 


DlSiZzJ, 6144,16 


! 12 SCANS OF 1024 8 2*8 BYTES PER RECORD 


270 


CREATE BDAT 


D2$6Zz$,6144,16 


! 1 024*1 2/2=NUHBER OF REC0RDS=6144) 


280 


PRINT USING 


'V 




299 


PRINT * 


DATA WILL BE STORED ON DISKETTE WITH FOLLOWING FILE NAMES: 1 


300 


PRINT USING 


'll/' 




319 


PRINT * 




NO PARTICLES — - FILENAME = ";D1$ 


320 


PRINT D 




PARTICLES FILENAME = ";D2$ 


339 


PRINT USING 


'll' 




340 


PRINT * 


IS THIS A CALIBRATION ? ENTER 7 Y 7 IF YES 1 


350 


PRINT “ 




ANYTHING ELSE IF NO 1 


360 


INPUT R$ 






370 


PRINT USING 


'//' 




380 


PRINT ' 




BE SURE LASER IS ON* 


390 


PRINT " 




PRESS [CONTINUE] WHEN READY* 


400 


PAUSE 
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410 I c SS=‘Y’ THEN 470 
420 PRINT USING 

430 INPUT 'ENTER THE THRESHOLD PRESSURE TO TRIGGER THE DEVICES (psi)\S8 

440 INPUT 'ENTER TINE DELAY FROM THRESHOLD PRESSURE (sec) ',T8 

450 ! TIMES ARE INTERPRETED BY THE COMPUTER IN SECONDS DOWN TO .001 

460 V0=S3/151.5 ! CALIBRATION OF PRESSURE TRANSDUCER 

470 D3$=D 1 $ ! NO-PARTICLES STRING NAME 

480 GOSUB Hultiprog 

490 GOSUB Storedata 

500 IF R$=*r THEN 1000 

510 |**x*******x********************************x**************************** 

520 Contcheck: IContinuity checking 

530 OUTPUT 709;*AC20‘ ICONNECT 3456A DVM TO IGNITION FIRING LINE 

540 LOCAL LOCKOUT 7 'KEEPS OPR. FROM SWITCHING FRONT PANEL TO LOCAL 

550 R 7=225 ! IN THOUSANDS OF OHMS 

560 OUTPUT 722; 'HSM002SW2SO1L1S0F4R1T35X1 * IRESISTANCE MEASUREMENT 

570 ! H RESETS DVM;SM002 SETS SERVICE REQUEST MASK WHERE 002 IS OCTAL REP. OF 

580 ! THE SERIAL POLL MASK BYTE;SW2 TELLS WHICH TERMINAL SWITCH IS USEDjSOl 

590 ! SYSTEM OUTPUT MODE ON - WAITS FOR CONTROLLER TO HANDSHAKE ;L1 LOAD 

600 ! INTERNAL MEMORY ON;SO-FUNCTION SHIFT OFF;F4-TWO WIRE CONNECTION TO DVM 
610 ! R1 -AUTORANGING ;T3-SINGLE TRIGGER ;Q-LOAD INTERNAL MEMORY OFF jXl -EXECUTE 
620 ! STORED PROGRAM, 

630 GOSUB Reading 

640 RS=V 

650 IF R8<R7 THEN GOTO Ok 

660 |******x******x********CONTINUITY CHECK********************************* 

670 PRINT USING '8' 

680 PRINT 'CONTINUITY CHECK BAD!!!!' 

690 PRINT 'RECHECK BEFORE PROCEEDING. WHEN CHECKED, PRESS [CONTINUE! ' 

708 BEEP 3000, .3 

710 WAIT .1 

720 BEEP 100,1.0 

730 PAUSE 

740 GOTO Contcheck 

750 OkiPRINT USING '8' 



760 


DISP ' 


CONTINUITY CHECK IS O.K.' 


770 


PRINT ' 


BE SURE NITROGEN IS ON' 


780 


OUTPUT 709 j 'AC21 ' 


ICONNECT DVM FOR IGNITOR VOLTAGE MEASUREMENT 


790 


PRINT ' 


DVH CONNECTED TO FIRE SWITCH' 


800 


PRINT USING V 




810 


PRINT 'BE SURE VISICORDER 


IS SET UP TO RUN ON PROPER SCALE WITH LAMP ON.' 


CO 

ro 


PRINT USING '/' 




830 


OUTPUT 722 ; 'HSM0 02SU2Z0SO1 L1FL0 .01STIS0F 1R1 T3Q ' ! VOLTAGE ON IGNITER 


840 


DISP ' 


STANDING BY FOR IGNITION' 


850 


PRINT * 


STANDING BY FOR IGNITION* 


860 


BEEP 2000,. 1 




870 


OUTPUT 722; 'XI' 


! VOLTMETER TRIGGER 


880 


GOSUB Reading 


.'READS VOLTMETER 


890 


R9=ABS(V) 




900 


IF R9< 1 0 THEN GOTO 870 


! 12 VOLTS ON IGNITOR 
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910 O0=TIMEDATE 
920 OUTPUT 709; "AC22 ' 

930 OUTPUT 722; *Xr 

940 GOSUB Reading 
959 R9=ABS(V) 

961 IF R9< VO THEN GOTO 930 
970 WAIT T8 
980 01=TIMEDATE 




! TRIGGER VOLTMETER 



! THRESHOLD PRESSURE 
! TINE DELAY 



998 IF R$( )*Y‘ THEN 1060 
1 000 ! 



! IF ACTUAL RUN THEN SKIP SOME LINES 
FOR CALIBRATIONS INTRODUCE PARTICLES AND THEN CONTINUE 



1010 PRINT CHR$(12) 

1028 PRINT USING ’////////’ 

1030 PRINT * INTRODUCE PARTICLES THEN PRESS KEY 4 9 TO TAKE DATA’ 

1040 ON KEY 9 LABEL * TAKE DATA* GOTO 1060 
1050 Standby: GOTO Standby 

1060 GOSUB Hultiprog ! TAKE PARTICLE DATA 

1070 LOCAL 7 

1080 D3$=D2$ ! PARTICLE DATA FILE NAME 

1090 GOSUB Storedata 

1180 PRINT * ELAPSED TIHE: FIRE TO MULTIPROGRAHHER CALL = * ‘,01-00 

1110 PRINT * ELAPSED TIHE: FIRE TO MEMORY INTERRUPT = *:02-00 

1120 PRINT * DATA STORED ON DISK WITH FILENAMES C;D1$;*) AND ( , ;D2$; , ) > 

1130 GOTO End 

1149 Hultiprog; ! 

1150 OUTPUT 723; “CC, 2, 3,12, 13T* 

1160 OUTPUT 723j"CC,5,6,9 f 10T* '.CLEARS THE ARM, BUSY AND EOP OF MEMORY CARDS 
1170 OUTPUT 723; 'CC, 1 ,4, 11T* .'CLEARS SAME FOR A TO D'S AND TIMER PACER 

1180 OUTPUT 723;*SF, 2, 3,1,. 001,12, T* ! THE (1) IS 2'S COMPLIMENT BINARY 

1190 OUTPUT 723; *SF,5,3, 1 , .001 , 12,T* ! THE .001 IS THE LEAST SIGNIFICANT BIT 

1200 OUTPUT 723;*SF,9,3,1 , .001, 12, T* ! THE 12 IS FOR 12 BIT WORD SIZE 

1210 OUTPUT 723; *SF, 12, 3,1, .001, 12, T* ! SINCE THE A TO D IS 12 BIT 

1220 OUTPUT 723 ; ’WF , 3 , 1 023T , WF , 6 , 4 095T ’ 

1230 OUTPUT 723j , WF,1G,4095T,UF,13,409jT* 

1240 ! SETS REFERENCE WORD FOR WHEN TO STOP TAKING DATA AND GENERATE INTERRUPT 
1250 OUTPUT 723 ; *UF , 2 . 1 , 1 T , UF , 5 . 1 , 1 T , UF , 9 . 1 , 1 T , WF , 1 2 . t , 1 T * 1SETS FIFO MODE 
1260 OUTPUT 723;*WF,3, 1 ,0T,WF,3.2,0T,WF,3.3,0T,WF,13. l,0T,WF,13,2,OT,WF,13.3,0T 



1300 OUTPUT 723; *UF,4 ,2,0T,UF,4,1ST- ! TIMER PACER GIVES 1 PULSE OF 1 SEC WHEN 
1310 ! TRIGGERED BY THE BLANKING PULSE (PLENTY OF TIME FOR 8 SCANS) 
1320 WAIT 2 

1330 K=SP0LL(723) ! WAIT FOR HEM INTERRUPT 
1340 IF K(>64 THEN GOTO 1330 

1350 OUTPUT 723; *WF ,4.2, 1T“ ! MAY BE UNNECESSARY TO ALTER TIMER PACER SINCE 
1360 ! MEMORY CARDS HAVE AUTOMATIC LOCKOUT BUT FOR NOU WE WILL DO IT 
1370 02=TIMEDATE ‘ 

1380 ON ERROR GOTO Err trap ! NEEDED TO READ ARMED CARD INTERRUPT LIST 
1390 SEND 7 ; UHL MLA TALK 23 SEC 12 ! SPECIFICALLY ASKS FOR INTERRUPT LIST 




1280 



! SETS COUNTER AND POINTERS OF 2ND MEMORY CARD IN EACH PAIR TO 0 
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1430 y a r _ r e a c : ENTER 7 , Addr 9 

1410 fleficaros: PRINT ’ 

1420 PRINT 1 

1420 HAT Address= (0) 



3 (») ! READ WHICH CARDS INTERRUPTED 

MEMORY CARDS WHICH GENERATED INTERRUPTS ARE 1 
SLOTS* = * j Address ( * ) 



1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1599 

1600 
1610 
1620 



OFF ERROR 

OUTPUT 723j*DC,3,6,13 > 10T* 
OUTPUT 723; “HI,2,1024T “ 
ENTER 72305; A( *) 

OUTPUT 723;'HI,5,4096T” 
ENTER 72305;B(») 

PRINT * 

OUTPUT 723; , HI,12,4096T' 
ENTER 72305;C( *) 

OUTPUT 723;*HI,9,4096T* 
ENTER 72305;D(») 

ENABLE INTR 7;8 
PRINT 1 
HAT A= (-1 )*A 
HAT B= (-1XB 
HAT C= (-1 )*C 
HAT D= (-1 )*D 
RETURN 
Storedata; ! 



EXHAUST 



1DISSARH MEM CARDS 
! SET UP CARD TO BE READ 
GETS DATA FROH 1024 HEHORY BOARD 



! GETS DATA FROH 4096 HEHORY BOARD 
DATA ENTERED* 



1 



MOTOR DATA ENTERED* 

THE 1024 CARD IS INCLUDED BUT NOT SAVED, IT DIDN'T 
PERFORH WELL, COULD BE REPLACED BY A 4096 CARD. 



! DIODE VOLTAGES ARE NEGATIVE SO SIGNS ARE CHANGED 



1630 ASSIGN PDiskf ile TO D3HZz$ 

1640 ASSIGN g Buffi TO BUFFER C(«) 

1650 ASSIGN ?Buff2 TO BUFFER DU) 

1660 ASSIGN gBuff3 TO BUFFER B(*) 

1670 CONTROL §Bt»ff 1,3; 1,32768,1 
1680 CONTROL §Buff2,3;l, 32767,1 
1690 CONTROL 8Buff3,3;l,32767,l 
1700 TRANSFER SBuff 1 TO gDiskfile 
1710 WAIT FOR EOT PDiskfile 
1720 TRANSFER gBuff2 TO gDiskfile 
1730 WAIT FOR EOT gDiskfile 
1740 TRANSFER gBuff3 TO gDiskfile 
1750 WAIT FOR EOT gDiskfile 
1760 ASSIGN gDiskfile TO * 

1770 ASSIGN gBuffl TO * 

1780 ASSIGN gBuff 2 TO * 

1790 ASSIGN gBoff3 TO * 

1800 RETURN 
1810 Reading : ! 

1820 STATUS 7,1 ;A8 
1830 ENTER 722;V 
1840 ENABLE INTR 7;8 
1850 RETURN 

1860 Errtrap : IF ERRN=159 AND ERRL(Var 

1870 PRINT ERRH$ 

1880 GOTO Hencards 
1890 End: END 



! HOTOR CAVITY 4 SCANS 
! HOTOR CAVITY 4 SCANS 
! EXHAUST 4 SCANS 
! SETS BUFFER POINTERS TO FULL 
'INTERFACE REGISTERS SECTION OF 
'LANGUAGE HANUAL 
! ORDER OF DATA ON THE DISX IS 
! HOTOR CAVITY — 8 SCANS 
! EXHAUST -- 4 SCANS 
'WAIT BECAUSE OVERLAPPING 
'TRANSFERS ARE NOT WANTED 

'CLOSE I / O PATHS 



! CHECKING STATUS BEFORE READING 
•VOLTMETER IS A FORHALITY TO 
'CLEAR THE SERVICE REQUEST 

read) THEN Hencards 



'EVEN IF THE ERROR WAS NOT THE ONE PLANNED 
!FOR PROGRAH EXECUTION CONTINUES 
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10 i ******************************************************************* 



20 !******************* 

30 j******************* 

40 i******************* 

50 j******************* 

60 !******************* 

70 !******************* 

80 !******************* 

90 !******************* 

100 I******************* 

118 
120 
130 
140 
150 
160 
R 

170 
180 
198 
200 

210 Choose: 

228 01d=O 

230 PRINT 
ON’ 

240 PRINT 
250 PRINT ’ 

260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 



RDC 

PLOTS RAW DATA 
AVERAGES 
FILTERS 

DETERMINES MEAN DIAMETER 
BY INTERACTIVE GRAPHICS 
AND THE TWO-ANGLE METHOD 
Robert Kelly Harris 
1984 



******************** 

******************** 

******************** 

******************** 

******************** 

******************** 

******************** 

******************** 

******************** 



!******************************************************************* 

OPTION BASE 1 

COM /Hrdoaus/ Av2(1024) 

COM /Gauss/ T1U 024),G(1024) >L 
COM /Max/ M7,M5,Xt,Yt,X*,Yn 

COM /Readata/ B,P,H,Q3$l20],Q4Sl20],Zz*l20],Yl(8192) BUFFER, Y2( 81 92) BUFFE 

COM /Two/ Avl(1024),M,Ml,F 
DIH Scans(8),X( 1024) 

COM /Plots/ PHI201 ,P2$I201 ,P3$120I ,P4$1201 
INTEGER Graft 1:12480) BUFFER 
PRINT CHR S ( 1 2) 

* TO LOOK AT ANY DATA FILE THE PROGRAM NEEDS SOME STARTING INFORHATI 

• * 

TO ACCOUNT FOR THE CHANGE IN WAVELENGTH IN THE MEDIUM* 

PRINT ” 

PRINT * ENTER THE INDEX OF REFRACTION OF THE MEDIUM* 

PRINT " 

PRINT* WATER =1,33* 

PRINT * AIR = 1.8“ 

PRINT’ ESTIMATE OF EXHAUST = 1.1’ 

BEEP 

INPUT * THIS VALUE ADJUSTS THE COMPARISON CURVE TO THE MEDIUM’, Ml 

PRINT CHR$( 12) 

PRINT * TO ACCOUNT FOR REFRACTION OF LIGHT AT BOUNDARIES BETWEEN* 

PRINT ** 

PRINT * THE MEDIUM AND AIR YOU MUST ENTER THE INDEX OF REFRACTION * 
PRINT *“ 
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400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 



Cr THE COMBINATION Or THE MEDIUM AND ITS BOUNDARY' 
THIS VALUE IS APPLIED DIRECTLY TO THE DATA' 



PRINT ' 

PRINT " 

PRINT 1 
PRINT ” 

PRINT 1 AIR=1 ,0 PLEXIGLASS 6 WATER = 1.39' 

PRINT ’* 

PRINT 'ESTIMATE OF MOTOR CAVITY COMBUSTION PRODUCTS 6 WINDOW = 1.22' 
PRINT ,B 

PRINT ' ESTIMATE OF ROCKET EXHAUST = 1.1 OR 1,0' 

INPUT M 

PRINT CHR$( 12) 

INPUT * ENTER LASER WAVELENGTH (HeNe=, 6328, Ar=. 488)', L 

INPUT 'ENTER THE DESIRED PLOTTING INTERVAL OF DIODE ARRAY (2,4,6)\H 



F=460 
Diod=.025 

FOR 1=0 TO 1023 
X(I+1)=I+1 

Tl(I+l)=(I*Diod)/(F*M) 

NEXT I 

FRINT CHR$(12) 

CALL Displavl (Old) 

IF Old THEN CALL Review (Av 1 (*) ) 

IF Old THEN GOTO Gauss 
BEP 

INPUT ' INPUT TWO FILENAMES NOW (NO-PART ,PART) ',B3$,Q4$ 

INPUT ' ENTER 'V FOR MOTOR '2' FOR EXHAUST', P 

Zzt=' : INTERNAL ,4,1“ ! USES LEFT DISK DRIVE FOR RAW DATA FILES 

CALL Readata 

CALL Shift 

BEEP 



!hm FOCUL LENGTH OF OBJECTIVE LENS 
!n« DIODE SPACING HAY BE .03 

(CREATE AN ARRAY OF DIODE NUMBERS 
(CREATE AN ARRAY OF THETA 



(JUST PRINTS INITIAL REMARKS ON CRT 
(THESE TWO LINES ALLOW FOR 
(REVIEWING REDUCED DATA FILES 



690 Screen: ! 

700 PRINT CHRSU2) 

710 PI $= 'VOLTAGE * 

720 P25='DI0DE NUMBER' 

730 P3i=' VOLTAGE vs. DIODE* 

740 P4$=’RAW DATA' 

750 H5=1024 

760 Xt=64 

770 Yt=.04 

780 Xn=4 

790 Yn=5 

800 CALL Plot(Y2(») ) 

810 GSTORE Graft*) 

820 CALL Dataplot (B,Y1 (*) ,H) 
830 PRINT ' 

840 PRINT ” 

850 PRINT ' 

860 INPUT J 



(CLEARS SCREEN 
! STRINGS FOR PLOTS 



! SET UP VALUES FOR PLOTS 



! DRAWS AXES ,ETC... 

! STORES GRAPHICS DISLAY JUST MADE 
(PLOT NO-PARTICLES DATA 
NO-PARTICLES PLOTTED' 

HOW MANY SCANS SEEM TO BE GOOD ?' 
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370 PRINT ■* 

385 PRINT * WHICH SCANS ARE GOOD ?. . . 10 , . ,1 ,2, 4, 5, 7, INCLUDE LAST CO 

MMA* HAST COMMA IS REQUIRED OR YOU HAVE TO HIT CONT,, TWICE 

890 INPUT Scans!*) 

900 PRINT CHRS! 12) 

910 PRINT USING •/////////■ 

920 PRINT 1 AVERAGING THE SELECTED SCANS' 

930 GRAPHICS OFF 

940 CALL Averaged, Scans!*) ,Y1(») ,Avl(*)) 

950 PRINT USING *g' 

968 GLOAD Graf 00 ’LOADS GRAPHICS ARRAY RATHER THAN WASTE TINE RE-DRAWING 
970 GRAPHICS ON 

988 CALL Result(Avl(*),X(*),H) ! PLOTS 1024 ELEMENT ARRAYS 

990 PRINT USING '//' 

1000 PRINT * Average Intensity No-Particles' 

1010 ON KEY 0 LABEL 'AVERAGE 1 ' GOTO Screen 
1020 ON KEY 1 LABEL 'PLOT-PARTICLES' GOTO 1060 

1030 PRINT * PRESS KEY * 0 TO RE-AVERAGE OR * 1 TO CONTINUE' 



1040 Standby: 



! MANUAL CALLS THIS INTERRUPT DRIVEN PROGRAMMING 
! LOOPS, WAITING FOR USER TO DECIDE 



'HELPS AVOID CONFUSION BY CLEARING THAT BOX 

! PLOT PARTICLES DATA 



1050 GOTO Standby 
1060 PRINT USING *?" 

1070 OFF KEY 1 

1080 CALL Dataplot(B,Y2(*),H) 

1090 PRINT USING 7* 

1100 PRINT ' PARTICLE DATA PLOTTED* 

1110 PRINT * FOR A HARD COPY OF THIS RAW DATA PRESS KEY 8 6' 

1120 PRINT * ' 

1130 PRINT * OR TO CONTINUE PRESS KEY * 1' 

1140 ON KEY 6 LABEL 'HARD 4 RAW* GOTO Raw 
1150 ON KEY 1 LABEL 'CONTINUE 1 GOTO Select 
1160 GOTO Standby 

1170 Raw: CALL Plot(Y2(*),l) ! THE ONE (1) IS AN OPTIONAL PARAMETER 

1180 CALL Dataplot(B,Y2(*),H) ! WHICH IS USED TO GET HARD COPIES 

1190 Select: ! 



-.1200 OFF KEY 6 

1210 OFF KEY 1 

1220 PRINT CHR$( 12) 

1230 PRINT USING '////' 

1240 PRINT * HOW MANY SCANS SEEM TO BE GOOD ?' 

1250 INPUT J 

1260 PRINT* WHICH SCANS ARE GOOD ?... ie... 1,2, 4, 5, 7, INCLUDE 

LAST COMMA* 

1270 INPUT Scans!*) 

1280 PRINT CHRS! 12) 

1290 PRINT USING */////////' 

1300 GRAPHICS OFF 

1310 PRINT * AVERAGING THE SELECTED SCANS' 

1320 CALL Average (J,Scans(*) ,Y2(*) ,Av200 ) 

1330 PRINT USING '?' 

1340 P 4$=’ AVER AGED SCANS' 
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1350 


CALL Plot(4v2<*)) 




1360 


GRAPHICS ON 




1370 


GSTORE Graft*) 




1380 


CALL Result(Av2(*),X(*),H) 




1390 


PRINT USING *//* 




1400 


PRINT * 


Averaae Scattered Intensity* 


1410 


ON KEY 2 LABEL 'AVERAGE2* GOTO 1080 


1420 


ON KEY 3 LABEL ‘SUBTRACT* GOTO 1450 


1430 


PRINT * 


HIT KEY 4 2 TO RE-AVERAGE OR KEY 4 3 TO 


CONTINUE* 




1440 


GOTO Standby 




1450 


PRINT USING *g* 




1460 


OFF KEY 3 




1470 


NAT Av 1= Av2-Av 1 


! SUBTRACTS NO-PARTICLES FRON PARTICLES 


1480 


GLOAD Graft*) 




1490 


CALL Result(Avl(*),X(*),H) 




1500 


PRINT USING *///* 




1510 


PRINT * 


Plot of the Difference Between Particles and N 


o-Particles* 




1520 


PRINT * 


KEY 4 6 FOR HARD COPY* 


1530 


PRINT * 


KEY 4 1 TO CONTINUE B 


1540 


ON KEY 6 LABEL * HARD AVERAGE* GOTO 1570 


1550 


ON KEY 1 LABEL * FILTER 


* GOTO 1600 


1560 


GOTO Standby 




1570 


CALL PlottAvl (*), 1 ) 




1580 


CALL Result(Avl(*),X(*),H) 




1590 


PLOTTER IS 3, ‘INTERNAL* 




1600 


PRINT USING *//* 




1610 


PRINT * 


ENTER THE NUNBER OF TINES YOU WISH TO APPLY* 


1620 


PRINT * 


THE DIGITAL FILTER FOR SMOOTHING* 


1630 


PRINT * 


EXAMPLE******** 10* 


1640 


PRINT * 


TAKES ABOUT 1.5 MINUTES* 


1650 


INPUT * 


ZERO IS ALSO ACCEPTABLE* ,Fil 


1660 


IF Fil=0 THEN Gauss 




1670 


P4$=*FILTERED DATA* 




1680 


CALL Filter(Avl(*),Fil) 




1690 


PRINT CHR $(12) 




1700 


CALL Plot (Avl Ot) ) 




1718 


IF Old THEN NAT X= T1 




1720 


CALL Res»lt(Avl(*),X<*),H) 




1730 


PRINT * 


Plot of the Difference Between Particles and N 


o-Particles* 




1740 


PRINT * 


AFTER APPLICATION OF A DIGITAL FIL 


TER* 






1750 


PRINT * 


for a HARD COPY * 


1760 


PRINT * 


PREPARE the PLOTTER and* 


1770 


PRINT * 


PRESS KEY 4 8* 


1780 


PRINT * 


OR* 


1790 


PRINT * 


PRESS KEY 4 9 TO NORNALIZE* 


1800 


ON KEY 8 LABEL 'HARD FILTERED* GOTO 1830 


1810 


ON KEY 9 LABEL * NORNALIZE 


' GOTO Gauss 
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1 The 1 is for hard copy 



1829 

1830 
1840 
1850 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2920 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2208 
2210 
2220 
2230 
2248 
2250 
2260 
2270 
2280 
2290 
2300 
2310 



GOTO S f snob v 
CALL Plot(AvlU) .1) 

CALL Resvlt(Avl!*),XU),H) 
OFF KEY 8 
Gauss: OFF KEY 2 
OFF KEY 8 

ON KEY 1 LABEL 1 FILTER 



!KEY 1 OPTION IS STILL VALID 



GOTO 1600 ! ALLOWS ONE TO SMOOTH OLD DATA 
! WHEN REVIEWING REDUCED FILES 

PRINT CHR$(12) 

!Av2 Array is Normalized BUT A v 1 ( * ) is Still Saved For Re-work if Needed 
MAT Av2= Avl ! TO BEGIN WITH THE ARRAY IS ASSUMED TO BE NORMALIZED 
! AND THE USER CHANGES THIS WITH INTERACTIVE GRAPHICS 
! IN THE SUBROUTINE ‘ConDore 1 
P 1 $= ’NORMALIZED INTENSITY 1 

P2$=’THETA (rad) 1 'STRINGS FOR PLOTS 

P3i=" INTENSI TY vs. THETA" 

P4t=' CURVE FIT RESULTS 1 
H5= . 05 
Xt= ,0 04 
Yt=.l 
Xm=2 
Yn=2 

CALL Plot < Av2( *)) 

CALL Result(Av2U),TlU),H> 

GSTORE Graft*) 

CALL Compare! Avl (*),H,M, Ml, Graft*)) 

GRAPHICS OFF 

ON KEY 4 LABEL ’OTHER ARRAY 1 GOTO Choose 
TWO-ANGLE" GOTO Bucbele 
NORMALIZE 1 GOTO Gauss 



! SET UP VALUES 
! FOR PLOTS 



! PLOT OF NORMALIZED INTENSITY vs. THETA 
! SAVES THE SCREEN IMAGE 



ON KEY 5 LABEL 
ON KEY 9 LABEL 



! PRINTS OPTIONS ON THE SCREEN 



ON KEY 6 LABEL 'STORE DATA 1 GOTO Storedata 
ON KEY 7 LABEL 1 QUIT 1 GOTO Quit 
CALL Menuel 
GOTO Standby 
Buchele: ! 

CALL Twoangle 
CALL Menuel 
GOTO Standby 
Storedata : ! 

CALL Store 
CALL Menuel 
GOTO Standby 
Quit: ! 

END 

SUB Averaoe!J, Scans!*), YU) ,Av(*)) ! AVERAGES SELECTED SCANS 

MAT Av= !3) UNITIALIZES THE ARRAY LOCAL TO THIS ROUTINE 

FOR 1=1 TO J STEP 1 

K=!Scans(I)-l)*1024+l ' ! THIS COUNTER IS THE BEGINNING 

L=Scans(I)*1024 ! AND THIS ONE THE END OF BLOCKS 

FOR I1=K TO L ! OF 1024 IN THE OVERALL DATA 
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2320 I2=IHH1 

2330 Av(I2)=4v(I2)+Y(I!) 

2340 NEXT I! 

2350 NEXT I 

2360 MAT Av= Aw/CJ) 

2378 SUBEND 

2380 SUB Dataplot(B,Y(*),H) 

2390 COM /Max/ M7,M5,Xt, Yt,Xn,Yn 

2480 LDIR 0 

2418 L$G 4 

2420 FOR J=1 TO B 



2480 MOVE 1,Y(K) 

2490 FOR I=K TO L STEP H 

2500 X=I-K+1 

2510 PLOT X, Y(I) 

2520 NEXT I 

2530 NEXT J 

2540 PENUP 

2550 SUBEND 

2560 SUB Result(Iav(*),D(*),H) 

2570 CLIP ON 

2580 MOVE D(l),Iav(l) 

2590 FOR 1=1 TO 1024 STEP H 

2600 DRAW D(I),Iav(I) 

2610 NEXT I 

2620 PENUP 

2638 SUBEND 

2648 SUB Plot (Y(») OPTIONAL Hard) 
2650 COM /Max/ M7,H5,Xt ,Yt,X«,Y« 

2660 COM 

2670 M7=MAX(Y(i))«l.l 
2680 H6=.1*M7 
2690 PRINT ” 

2780 GINIT 
2718 

2720 SELECT NPAR 
2730 CASE 1 
2740 PLOTTER 
2750 CASE 2 
2760 PLOTTER 
2770 END SaECT 
2780 GRAPHICS ON 
2790 CSIZE 5, .6 
2800 DEG 
2810 LDIR 0 



! THIS COUNTER 13 ALWAYS BETWEEN 
! 1 AND 1024 



!*****THE AVERAGE*****! 

!B IS 4 OR 8 (THE NUMBER OF SCANS) 



(EACH SCAN 
! SEE NEXT LINE 

KEEP TRACK OF EACH SCAN AS IT APPEARS 

(BEGINNING OF EACH SCAN AND 
! THE END WITHIN THE TOTAL BLOCK 
! MOVE TO THE FIRST POINT 

! THIS GIVES 1 TO 1824 FOR ABSCISSA 



! THE AVERAGE INTENSITY IS PLOTTED 



! IF THE OPTIONAL(Hard)IS RECEIVED 
•THE FIGURE GOES TO THE PLOTTER 
/Plots/ Pm20],P2M20],P3$l20],P4$[20] 

! A SCALING VARIABLE 
! ANOTHER SCALING VARIABLE 

! GINIT IS JUST GOOD PRACTICE SO YOU KNOW 
•WHERE YOU ARE BEGINNING 
!THIS DETECTS IF THE HARD COPY IS DESIRED 
• IT COULD BE DONE WITH IF THEN LOGIC 
IS 3, "INTERNAL’ !BUT IS PRESENTED FOR FAMILIARIZATION 

! SINCE IT IS MORE POWERFUL IN COMPLEX 
IS 705,’HPGL" ! SITUATIONS 

! TO BE ABLE TO SEE THE PLOT 



2430 MOVE 500+J*50,M7- .05 

2448 LABEL J; ! HELPS 

2450 LINE TYPE 1 

2460 K=(J-1)*1024+1 

2470 L=J*1024 
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2820 LORG 5 

2935 MOVE 75,95 

2840 LABEL USING ‘K‘;P4$ 

2950 HOVE 75,90 
2860 LABEL P3$ 

2870 HOVE 75,20 
2880 LABEL P2$ 

2898 LDIR 90 
2900 HOVE 18,60 
2910 LABEL Pl$ 

2920 VIEWPORT 30,125,30,85 
2938 FRAME 

2940 WINDOW 0,H5,-H6,H7 

2950 AXES Xt,H6,0,0,X*, 1,5 

2960 LDIR 0 

2970 CSIZE 3,, 5 

2980 LORG 8 

2998 CLIP OFF 

3000 FOR I=-M6 TO H7 STEP H7/10 

3010 HOVE 0,1 

3020 LABEL USING **,HD.DDDD*;I 

3030 NEXT I 

3040 CSIZE 3, .6 
3050 LDIR 90 
3060 
3070 
3088 
3090 
3109 
3118 



! SCREEN UNITS FOR MARGINS 
! DRAWS A MX 
! EXTENT OF X AND Y 



! NUMBER THE Y AXIS 



'NUMBER THE X AXIS 



LORG 8 

FOR 1=0 TO M5 STEP Xt 
MOVE I,-H6 
LABEL USING **,K’;I 
NEXT I 
PENUP 
3120 SUBEND 

3130 SUB C<mpare(Avl(*),H,H, Ml, INTEGER Graf(*)) 

3140 COM /Gauss/ Tl(»), GOO, L 

3150 COM /Hrdgaus/ Av200 

3160 COM /Max/ M7,M5,Xt,Yt,X«,.Yrt 

3170 DIM TO :1024),At(l : 17) ,Dobbins(l : 17) ,Tld(l :17) 

3180 DATA 3., 3. 5, 4., 4.5, 5., 5.5, 6., 6. 5, 7. ,7, 5, 8., 8. 5,9., 9, 5,18, ,10.5,11. 

3190 DATA .0556,, 035,, 0206,. 014,. 0106,. 0081,. 00605,. 005,. 00374,. 0032,, 00248, 
.0022, .00185, .0016, .00135, .0012, .001 
3200 READ AtOO 

3210 READ DobbinsU) 

3228 OFF KEY ! GETS RID OF ALL LABELS ON KEYS 

3238 ON KEY 3 LABEL ‘MENUE’ GOTO Subexit 

3240 PLOTTER IS 3, ‘INTERNAL* ! IN CASE A HARD COPY WAS JUST MADE 

3250 Hard=0 ! (0) SO ONE DOESN'T EXIT TOO SOON 

3260 Centerline=l ! THE INITIAL NORMALIZING VALUE 

3270 D=2Q . ! INITIAL PATICLE MEAN DIAMETER IN MICRONS 



3280 Change: CSIZE 4,. 6 

3290 PRINT USING 7 /////* 
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3300 

3310 

3320 

3330 

3340 

3350 

3360 

3370 

3380 

3390 

3400 

VE 

3410 

3420 

3430 

3440 

3458 

3460 

3470 



PRINT 1 
MOVE 0,1 
CLIP ON 
Con=PI«D*Ml/L 
FOR 1=1 TO 1024 STEP H 
T ( I )=T 1 ( I ) *C on 
G(I)=EXP(-((,57*T(I)) A 2)) 
IF T ( I ) ) 3 . 0 THEN 3400 
DRAW T 1 < I) ,G(I) 

NEXT I 

MAT Tld= At/(Con) 

FOR 1=1 TO 17 
DRAW Tld(I) ,Dobbins< I) 
NEXT I 
PENUP 

IF Hard THEN SUBEXIT 
PRINT USING •millin' 

PRINT * 



STAND BY FOR CURVE* 



! A CONSTANT 

! Theta-bar for sone Dbarlnote Ml) 
.'Gaussian curve for this Dbar 
! GAUSSIAN NO GOOD FOR Theta-bar )3 
.'Gaussian curve vs, theta 

ITHETA FOR POINTS ON DOBBINS'CUR 



! EXITS ROUTINE IF A HARD COPY WAS JUST MADE 
USE THE KNOB TO VARY THE PARTICLE 



SIZE 1 

3480 


PRINT ’ 


OR HIT KEY * 6 FOR HARD COPY’ 


3490 


PRINT ’ * 




3500 


PRINT ’ 


KEY * 9 ALLOWS YOU TO* 


3510 


PRINT " 


RE— NORMALIZE" 


3528 


PRINT ’ * 




3530 


PRINT ’ 


HIT KEY * 3 TO GET OUT* 


3540 


PRINT USING •l/ll/r 




3550 


ON KEY 6 LABEL ’HARD COPY 1 


GOTO Hardgauss 


3560 


ON KEY 9 LABEL ■NORMALIZE" 


GOTO Divide 


3578 


ON KNOB ,05 GOTO Pulse 


! ( , 05) IS TIME INTERVAL IN WHICH 


3580 Wait; GOTO Wait 


! PULSES FROM THE KNOB ARE 


3590 Pulse: PRINT USING ’S’ 


! COUNTED AND THIS NUMBER IS 


3600 


Strngl$ =, Mean Diameter * 


! USED BY THE INTERRACTIVE 


3610 


Strna2$=* Microns' 


! GRAPHICS TO VARY THE PARTICLE 


3628 


Count=KNOBX 


! SIZE AND PLOT THE ASSOCIATED 


3630 


D=DROUND( D+Coun t/15 ,2) 


! GAUSSIAN APPROXIMATIONS OF 


3640 


GLOAD Graf(x) 


! SCATTERING PROFILES 


3650 


LDIR 8 




366 0 


LORG 8 




3670 


MOVE H5,.8*H7 




3680 


LABEL StrnglUVAL$(D)4Strnq 


2$ 


3690 


GOTO Change 




3700 


Hardgauss: ! 




3710 


PRINT ' 


PREPARE the PLOTTER' 


3720 


PRINT * 


PRESS CONTINUE fur* 


3730 


PRINT ' 


HARD COPY' 


3740 


PAUSE 




3750 


Hard=l 


! SO THAT SUBEXIT OCCURS AFTER HARD COPY 


3760 


CALL Plot(Av2(*),l) 


! 1**HARD COPY 


3770 


CALL Result(Av2(*),Tl(*),H) 
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3780 


PRINT USING *?■ 




3790 


LORG 8 




3808 


HOVE H5, ,8*H7 




3810 


LDIR 0 




3820 


LABEL StrnqlHVALi(D)&Strng2 


$ 


3830 


HOVE 0,1 




3840 


GOTO Chanqe 




3850 


Divide: ! 




3880 


PEN -1 


(ERASES LAST PLOT 


3870 


HOVE 0,0 




3880 


FOR 1=1 TO 1024 STEP H 




3890 


DRAW T1 (I) ,Av2(I) 




3900 


NEXT I 




3910 


PRINT CHR$(12) 




3920 


PRINT USING ■//////■ 




3930 


PRINT ’ 


USE THE KNOB TO RE-NORHALIZE* 


3940 


PRINT 1 


KEY * 4 WILL RE-PLOT USING" 


3950 


PRINT ' 


UPDATED SCALE* 


3980 


ON KEY 4 LABEL 1 RE-PLOT* 


GOTO Replot 


3970 


ON KNOB ,05 GOTO Vary 




3980 


GOTO Wait 




3990 


Vary : PRINT CHR$(12) 




4000 


PEN -1 


(ERASES LAST PLOT 


4010 


HOVE 0,0 




4020 


FOR 1=1 TO 1024 STEP H 




4030 


DRAW T1 (I),Av2(I) 




4049 


NEXT I 




4050 


PEN 1 




4080 


C2=KN0BX 




4070 


Center line=Cen ter line+C2/120 


(KNOB USED TO VARY NORHALIZING VALUE 


4080 


HAT Av2= Avl/(Centerline) 


(NEW NORHALIZED DATA ARRAY 


4090 


GOTO 3870 


(PLOTS THE NEW NORHALIZED ARRAY 


4100 


Replot: ! 




4118 


PRINT CHR $(12) 




4120 


CALL Plot (Av2(* ) ) 




4130 


CALL Result(Av2(*),Tl(*),H) 




4140 


GSTORE Graft * ) 




4150 


GOTO Change 




4160 


Suberit: ! 




4170 


SUBEND 




4180 


SUB Henoel 




4190 


PRINT USING *?* 




4200 


PRINT * 


YOU CAN RE-NORHALIZE* 


4210 


PRINT ’*********************************-m*****-m*x***3***i*****x 


4220 


PRINT * 


YOU CAN RE-AVERAGE (New Data Only)* 


4230 


PRINT ’********************************x*******xm*x*m**xmxixx 


4240 


PRINT * 


PRESS KEY O 4 TO LOOK AT OTHER DATA* 


4250 


PRINT * 


OLD OR NEU * 


4269 


PRINT * 


MOTOR / EXHAUST* 


4270 


PRINT ’****************************X***X***X*#*#***#X**X*XX*I**X**X 
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4290 

4290 

4300 

4310 

4320 

4330 

4340 

4350 

4360 

4370 

4380 

4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4538 

4548 

4558 

4560 

4578 

4580 

4590 



PRINT 1 PRESS KEY * 5 ’ 

PRINT ’ FOR THE ' TWO- ANGLE ' METHOD" 

PRINT * OF PARTICLE SIZING’ ' 

PRINT ’XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXHXXiXXXXXHXXXXX’ 

PRINT ’ TO STORE THE REDUCED DATA PRESS KEY * 6 

PRINT , *XXXX»XXXXXXXXXXXXX*XX*XXXXXX**XX**XXXXXX*XH**XJ**X*XXXXX*** 

PRINT ’ TO END THIS SESSION PRESS KEY * 7’ 

SUBEND 

SUB Displayl (Old) 

PRINT ’ TO REDUCE NEW DATA PRESS KEY 4 1’ 

PRINT ” 

PRINT ’ TO REVIEW PREVIOUSLY REDUCED DATA PRESS KEY 4 2’ 

PRINT * ’ 

ON KEY 1 LABEL ’ NEW DATA’ GOTO New 
ON KEY 2 LABEL " OLD DATA’ GOTO Review 
Wait: GOTO Wait 
New: PRINT CHR$(12) 

PUT THE DISK IN THE LEFT DRIVE IF IT IS NOT ALREADY’ 
ENTER THE NAMES OF THE FIRST AND SECOND FILES,’ 

EACH FILE HAS DATA FROM BOTH DIODE ARRAYS,’ 

FIRST IS NO-PARTICLES 6 NEXT IS PARTICLES DH,D2$* 



! THIS VARIABLE IS PASSED TO THE MAIN PROGRAM TO INDICATE 
ITHAT THE DATA TO BE READ IS ONE (1) SCAN AND THAT 
‘.THE AVERAGING ROUTINES WILL NOT BE USED 

SUBEND 
SUB Readata 



PRINT ’ 

PRINT ” 
PRINT ’ 

PRINT ” 

PRINT ’ 

PRINT ” 

PRINT * 

PRINT ” 

Old=0 
SUBEXIT 
Review: 01d=l 



4600 

4610 

4628 

4630 

4640 

4650 

4660 

4670 

4680 

4690 

4700 

4710 

4728 

4730 

4740 

4750 

4760 

4770 



COH /Readata/ B,P.H.Q3$l2Q],Q4$l2Q],Zz$l20].Yl(#) BUFFER, Y2(x) BUFFER 



Xy z= 1 
B=8 

R 1=65536 

IF P=2 THEN Xy z=4097 
IF P=2 THEN R 1 =32768 
IF P=2 THEN B=4 
PRINT ’ " 

PRINT ’ READING 



! FILE POINTER VARIABLE 
! NUMBER OF SCANS IN MOTOR DATA 
INUMBER OF BYTES OF MOTOR DATA 
! RECORD 4 WHERE EXHAUST DATA BEGINS 
INUMBER OF BYTES OF EXHAUST DATA 
I NUMBER OF SCANS IN EXHAUST DATA 

DATA FROM FILE ON DISK’ 



ASSIGN gDiskl TO 93$6Zz$ 

ASSIGN gDisk2 TO Q4$6Zz$ IOPEN I/O PATHS 

ASSIGN gBuffl TO BUFFER Yl(x) 

ASSIGN gBuff2 TO BUFFER Y2(x) 

CONTROL gDiskl,5;Xyz 

TRANSFER gDiskl TO gBuffl ;COUNT R1 (READS NO-PARTICLE DATA 

WAIT FOR EOT gDiskl I NOT AN OVERLAPPING TRANSFER 

CONTROL gDisk2,5;Xyz IXyz SETS DISK FILE POINTER TO 

(EITHER MOTOR OR EXHAUST DATA 
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(READS PARTICLE DATA 



4783 TRANSFER «Disk2 TO §Buff2;C0UNT Rl 



4790 WAIT FOR EOT 9Disk2 
4990 ASSIGN iDiskl TO * 

4810 ASSIGN 8Disk2 TO * 

4820 ASSIGN gBuffl TO * 

4830 ASSIGN SBuff2 TO * 

4840 SUBEND 

4850 SUB Plot2(E,D(»),X> 

4860 DIN Titlel$[25),Theta$I20I 

4870 GINIT 

4880 IF X=1 THEN PLOTTER IS 705,*HPGL' 

4890 DEG 

4900 GRAPHICS ON 

4910 VIEWPORT 35,125,30,85 

4920 Hax=10*INT((HAX(D(*))+10)/10) 

4930 WINDOW D(1,2),D(E,2) ,0,Hax 

4940 G=(D(E,2)-D(1 ,2))/(E-l)*4 

4950 F=(INT(E/4)-l)*4 

4960 IF F=0 THEN F=4 

4970 CLIP D(4,2)-G,D(F,2)+2*G,0,Hax 

4980 GRID G,5,D(4,2)-G,9 

4990 CLIP OFF 

5000 LORG 8 

5010 LDIR 90 

5028 CSIZE 4, .5 

5030 FOR 1=4 TO E STEP 4 

5040 HOVE 0(1,23,0 

5050 LABEL USING '.DDDD'jDU, 2) 

5060 NEXT I 

5070 LDIR 0 

5080 LORG 8 

5090 FOR 1=10 TO Hax STEP 10 

5100 HOVE D(4,2)-G, I 

5110 LABEL USING , t,K , ;I 

5120 NEXT I 

5130 CSIZE 6,. 6 

5140 Titled 'TWO-ANGLE HETHOD* 

5150 Titlel$='For Various Angle Ratios' 
5160 Size*='SIZE (nicrons)' 

5170 Thta$='THETA (rad)* 

5180 Sub$='l" 

5190 LDIR 90 

5200 LORG 5 

5210 B=D(4,2)-G-(D(E,2)-D(1,2))/10 

5220 HOVE B,Nax/2 

5230 LABEL Size$ 

5240 LDIR 9 

5250 A=(D(E,2)+D(1 ,2))/2 

5260 HOVE A,-Hax/4 

5270 LABEL Thta$ • 



(JUST GOOD PRACTICE TO CLOSE 
(I/O PATHS 

(DEGREES for LABEL DIRECTION 

(AN X GRID LINE EVERY 4th POINT 
! MAKES GRID UNIFORH 

(PUTS NUHBERS ON X AXIS 

! NUMBER Y AXIS 

(STRINGS FOR PLOTS 
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5230 LOSS 3 

5230 MOVE A,-W4 

5380 LABEL Sub$ 

5319 LORG 5 



5320 

5330 

5340 

5350 

5360 

5370 

5380 

5390 

5400 

5410 

5420 

5430 

5440 

5450 

5400 

5470 

5480 

5490 

5580 

5510 

5520 

5530 

5540 

5550 

5560 

5570 

5580 

5590 

5600 

5610 



MOVE A,Hai*l ,15 
LABEL Titles 
CSIZE 4.5 
MOVE A,Haxil.05 
LABa Tit lei S 
PENUP 
SUBEND 

SUB Distribution(D(») ,Tratio,E) 
LINE TYPE 1 
CLIP ON 

Tratio$=VAL$ (Trat i o ) 

MOVE D(1,2),D(1,1) 

FOR 1=1 TO E 

DRAW D(I,2),D(I,1) 

NEXT I 
LORG 4 
LINE TYPE 1 
CSIZE 4,. 3 

MOVE D(I-1 ,2),D(I-l,l)+2 
LABEL Tra tio$ 

PENUP 

SUBEND 



! ANGLE RATIO 

•THIS SUBROUTINE PLOTS THE 
.'PARTICLE SIZES DERIVED USING 
(VARIOUS ANGLE RATIOS APPLIED 
! TO THE DATA OVER A RANGE OF 
(ANGLES 



SUB Twoangle 

( ***************************************************************** 



(*************** 
(*************** 
(*************** 
(*************** 
i *************** 



SUBPROGRAM 'TWOANGLE' 
PARTICLE SIZING BY 
MEASURING THE RATIO 
OF INTENSITY AT 
TWO ANGLES 



**********x******** 

******************* 

**********x******** 

******************* 

******************* 



(***************************************************************** 



5620 OPTION BASE 1 

5630 COM /Two/ Avl(1024),M,Hl,F 

5640 COM /Gauss/ T1(1024),G(1024),L (THE GAUSS IS NOT USED HERE 

5650 (BUT THE COM BLOCK HAS THETA 



5660 DIM D(200,2) (AND WAVELENGTH 

5670 PRINT CHR$(12) 

5680 PRINT USING *///////' 



5690 PRINT 
5708 PRINT 
5710 PRINT 
5720 PRINT 
5730 PRINT 
5740 PRINT 
5750 PRINT 
5760 PRINT 
5770 PRINT 



■ THIS SUBPROGRAM USES THE TWO-ANGLE METHOD DESCRIBED BY BUCHELE* 

I H 

* TO CALCULATE PATICLE SIZE FOR VARIOUS ANGLE RATIOS AND ANGLES.' 

I ■ 

1 IT IS HOPED THAT THE CURVES WHICH RESULT WILL SHOW 1 

■ • 

’ WHICH SIZE IS THE MOST PROBABLE.' 

■ • 

'AFTER NOTING FROM THE RAW DATA WHICH ANGLES CONTAIN THE CENTER LOBE' 
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5780 PRINT 
5790 PRINT 
5800 PRINT 
5810 PRINT 
5820 PRINT 
5830 PRINT 
5840 INPUT 



ENTER*************THE SMALLEST USEABLE ANGLE,* 

THE SMALLEST ANGLE RATIO, AND* 
THE STEP SIZE BETWEEN ANGLE RAIOS* 
YOU WISH TO EXPLORE* 
EXAMPLE************* .012, 1.2,. 4* 

ENTER thetal , theta-rat 10 , theta-step*,Q,A ,B 



5850 OFF KEY 

5860 X=0 ! GRAPH WILL APPEAR ON SCREEN 

5870 Begin: IRUNNING CONTINUES HERE WHEN A HARD COPY IS DESIRED 

5880 ! Van De Hulst and Gur.prech 6 Sleepevich explain that the chanae in 

5870 ! wavelength of the bean is accounted for by dividing by the index 

5700 lot refraction of the nediuri. *********L=L/M1 ********** 

5710 C=(L/Ml/.57/PI)*2 I see page 15 of nasa tech paper 2156 

5720 FOR N=1 TO 1024 ! t o see this is a convenient constant 



5730 IF 8<T1(N) THEN 5950 
5940 NEXT N 

5750 FOR Tra t i o=A TO 3 STEP B 



5960 
5970 
5980 
5990 
6000 
6018 
6020 
6030 
6040 
6050 
6060 
6070 
6080 
6090 
6100 
6110 
6120 
6138 
6140 
6150 
6160 
6170 
6188 

6190 Xcoup: 
6200 
6210 
6220 
6230 



FOR J=N TO 1 010/Tratio STEP 10 



'FINDS POSITION OF MINIMUM ANLGE 

! VAR IOUS ANGLE RATIOS 
! SET S THE RANGE OF POSSIBLE 
! ANGLES (THETA 1) 



Thta l=Tl (J) 

Thta2=Tl ( J*Tratio) 
Deltatheta=Thta2*2-Thtal *2 
11=0 
12=0 

FOR I=J-5 TO J+5 
I 1 =1 1+Av 1(1) 

NEXT I 

11 = 11/11 

FOR I=INT(J*Tratio-5) TO INT(J*Tratio+5) 



! THETA 2 



ITHESE 2 LOOPS AVERAGE A FEU 
! INTENSITIES IN THE HOPE 
! OF A STEADIER CALCULATION 



I2=I2+Avl(I) 
NEXT I 

12 = 12/11 
E=(J-N)/1 0+1 



! DONE HERE FOR THETA2 FOR 
! THE GIVEN ANGLE RATIO 



! THIS IS A COUNTER FOR THE ARRAY CONTAINING 
IPARTICLE SIZE AND THETA1 FOR THE GIVEN 
! ANGLE RATIO 

IF 11(12 OR 11=0 OR I2<=0 THEN D(E,1)=0 IALLOUANCE FOR IF THE DATA 
IF 11(12 OR 11=0 OR I2<=0 THEN GOTO Xcoup IIS NOT UELL BEHAVED 
D (E , 1 )=SOR (-C/Del tatheta*LOG( 12/11 )) ITUO ANGLE METHOD 

•VALUE OF DIAMETER based or. INTENSITY RATIO 
IFor a given ANGLE RATIO and ANGLE THETA1 

D(E,2)=Thtal 
Spar=PI*D(E,l )*M1/L 
Tbar=Spar*Thta2 
IF Tbar >3 THEN J=1 01 0/Tratio 
IF E=1 AND Tbar >3 THEN 6330 
6240 ! the above line ends all calculation if 

6250 ! the first elenent (E=l ) fa i led the test 
6260 NEXT J 

6270 ! FIRST TIME THROUGH- D HAS THE MOST ELEMENTS IT WILL HAVE 



! THIS IS PI*D/LAHBDA - THE SIZE PARAMETER 
ITHETA BAR FOR THE LARGER ANGLE 
'.THIS ENDS THE DO LOOP FOR THIS 
! ANGLE RATIO SINCE THE GAUSSIAN 
IIS NOT VALID WHEN Tbar ) 3 
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£0 SET UP THE GRAPH USING D's PRESENT PARAMETERS 



6230 ! 

6290 PRINT CHRSU2) 

6399 IF Tratio=A THEN CALL Plot2(E,D(*),X) 

6310 CALL Distribution(D(*),Tratio,E) 

6320 NEXT Tratio 



6330 ON KEY 2 LABEL ‘HARD COPY* GOTO Hard 

6340 ON KEY 3 LABEL ’HENUE * GOTO Sabexit 

6350 PRINT USING •/////■ 

6360 PRINT * YOU CAN GET A HARD COPY BY PRESSING KEY * 2‘ 

6370 PRINT * OR EXIT THIS ROUTINE BY PRESSING KEY t 3" 

6380 Standby :GOTO Standby 

6390 Hard :X=1 !A VARIABLE TO CONTROL PLOTTER 

6400 GOTO Begin 
6410 Subexit: GINIT 



6420 GRAPHICS OFF 

6430 SUBEND 
6440 SUB Shift 

6450 COM /Readata/ B,P ,H,Q3$[20I ,94$I20] ,Z:$I20J ,Y1 (*) BUFFER, Y2(») BUFFER 
6460 DIM E( 1:8192) 

6470 PRINT CHR$( 12) 

6480 PRINT USING "///////’ 

6490 PRINT * RAW DATA IS BEING SHIFTED TO CORRECT FOR SHALL GAPS BETWEEN SCANS’ 
6500 PRINT ,H 

6510 PRINT * BE PATIENT, IT'S A LONG SET OF DO-LOOPS* 

6520 ! THERE ARE SOHE SHALL GAPS BETWEEN SCANS AND THIS SUBROUTINE 

6530 ! SHIFTS THE DATA SO THAT THE FIRST DIODE DATA POINT IS HOVED 

6540 ! TO THE VERY BEGINNING OF ITS 1024 BLOCK IN THE OVERALL ARRAY 

6550 ! .THE FIRST SET IS RIGHT ON, THE NEXT IS ONE OFF, THE THIRD IS 

6560 ! TWO OFF ,,SO FORTH. THIS HAY NOT HATTER WITH OUR RESOLUTION 

6570 ! AND IS PROBABLY DUE TO THE HEHORY CARD CYCLING AT THE END OF 

6580 ! SOHE SCANS. SEE THE CIRCUIT TIHING DIAGRAM IN THE THESIS. 

6590 SELECT B 

6600 CASE 4 

6610 H=3 

6620 CASE 8 

6630 M=7 

6640 END SELECT 

6650 FOR K=1 TO 2 

6660 IF K=1 THEN HAT E= Y1 

6670 IF K=2 THEN HAT E= Y2 

6680 FOR J=0 TO H 

6690 IF H=7 AND J<=3 THEN 6800 !One 4096 Block Doesn't Need Shifting 

6700 FOR I=(J)*1024+1 TO (J+1 X1024 .'BLOCKS OF 1024 

6710 IF H=3 THEN L=I+J ! Array B Has the worst 

6720 IProblen With Shifting Data 

6730 IF H=7 THEN L=I+1 ! Array D is Allways off by one 

6740 IF LXJ+1 )*1024 THEN L=(J+1)*1024 

6750 ! JUST TO AVOID PROGRAMMING ERROR AT THE END OF THE ARRAY 

6760 E(I)=E(L) .'THIS SHIFTS THE PAT A 

6770 ' 'DEPENDING ON 'L', ARRIVED AT 



'EXHAUST DATA HAS 4 SCANS 



! MOTOR DATA HAS 8 SCANS 



'OPPERATES ON NO-PARTICLE AND PARTICLE SETS 
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!BY OBSERVING RAW DATA 



6780 

6790 NEXT I 

6800 NEXT J 

6810 IF K=1 THEN HAT Yl= E 

6820 IF K=2 THEN HAT Y2= E 

6830 NEXT K 
6840 SUBEND 

6850 SUB Store ! THE LARGE ARRAY CONTAINING HULTIPLE SCANS HAS BEEN REDUCED 

6860 ! TO A MEAN SCATTERING PROFILE BY AVERAGING AND FILTERING. IF 

6870 ! YOU FEEL THAT STORING THIS DATA IS NECESSARY THIS ROUTINE 

6880 ! DOES IT. IT SAVES DISK SPACE TO STORE THE RESULTS THEN 

6890 ! ELIHINATE THE RAW DATA IF YOU FEEL CONFIDENT THAT THE 

6900 ! REDUCTION IS THE BEST THAT CAN BE, IN OTHER WORDS, 

6910 ! DO NOT PURGE A RAW DATA FILE UNLESS YOU ARE ABSOLUTELY SURE 

6920 ! YOU WON'T WANT TO REDUCE IT AGAIN. 

6930 COH /Hrdaaus/ Av2(*) !Av2(* ) IS THE REDUCED DATA 

6940 DIH Data ( 1 : 1024) BUFFER ! BUFFER USED TO TRANSFER TO DISK 

6950 HAT Data= Av2 
6960 PRINT CHR$(12) 

6970 PRINT USING V/////* 

6980 PRINT * A SUGGESTED HETHOD OF NAMING REDUCED DATA FILES IS AS FOLLOWS* 

6990 PRINT * H HOTOR BEAM* 

7000 PRINT * X EXHAUST BEAM* 

7010 PRINT * C CALIBRATION , IF NO 'C' THEN AN ACTUAL FIRING IS ASSUMED* 

7020 PRINT * G GLASS BEADS * 

7030 PRINT * IF NO 'C' THEN 'G' STANDS FOR 'GAP' PROPELLANT* 

7040 PRINT * A ALUMINUM OXIDE* 

7050 PRINT * PPP CHAMBER PRESSURE FOR RUN OR OTHER. ... (NOZZLE TYPE)* 

7060 PRINT * PARTICLE SIZE FOR CALIBRATION* 

7070 PRINT * MMDD MONTH, DAY OF RUN OR CALIBRATION* 

7080 PRINT ** 

7090 PRINT * EXAMPLE: MCA125JN12* 

7100 PRINT ** 

7118 PRINT * HOTOR BEAH CALIBRATION USING ALUHINUH OXIDE FROH 1 TO 25 MICRONS* 
7120 PRINT " ON JUNE 12* 

7130- PRINT * EXAMPLE: XG550AU10* 

7140 PRINT ** 

7150 PRINT * EXHAUST BEAH DATA FOR GAP PROPELLANT AT 550 psi ON AUGUST 10* 

7160 PRINT * PLACE A DISK IN THE RIGHT HAND DRIVE* 

7170 PRINT * ENTER THE FILENAME YOU WISH TO USE FOR THIS REDUCED DATA* 

7180 INPUT At 

7198 CREATE BOAT At, 512, 16 

7200 ASSIGN BDisk TO At 

7210 ASSIGN BBuff TO BUFFER DataOO 

7220 CONTROL BBuff, 3;1, 8192,1 ! BUFFER IS FULL 

7230 TRANSFER BBuff TO BDisk ;COUNT 8192 
7240 WAIT FOR EOT gDisk 
7250 ASSIGN BBuff TO * 

7260 ASSIGN gDisk TO * 

7270 SUBEND 
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7280 

7290 

7300 

7310 

7320 

7330 

7340 

7350 

7360 

7370 

7380 

7390 

7400 

7410 

7420 

7430 

7440 

7450 

7460 

7470 

7480 

7490 

7508 

7510 

7528 

7530 

7540 

7550 

7560 

7570 

7580 

7598 

7600 

7618 

7620 

7638 

7640 

7658 

7660 

7670 

7680 

7690 

7700 



SUB ReviewtAvl (*)) 

DIM Datad : 1 024 ) BUFFER 

PRINT ‘ EACH REDUCED FILE CONTAINS ONE SCAN ONLY, YOU MUST KNOW IF IT 1 
PRINT * IS EXHAUST , MOTOR CAVITY , OR CALIBRATION DATA.* 

PRINT ” 

PRINT * THE DISK WITH THE REDUCED FILE SHOULD BE IN THE RIGHT-HAND DRIVE’ 
PRINT ” 

PRINT ’ ENTER THE FILENAME OF THE REDUCED DATA’ 

PRINT ’ TO BE REVIEWED’ 

PRINT ’’ 

INPUT A* 

ASSIGN gDisk TO At 

ASSIGN Muff TO BUFFER Data!*) 

CONTROL gDisk, 5;1 

CONTROL g B u f f , 3 ; 1 ,0,1 ITHIS IS AN EMPTY BUFFER 

TRANSFER gDisk TO Muff; COUNT 8192 !1024*8=NUMBER OF BYTES 

WAIT FOR EOT gDisk 
ASSIGN gDisk TO * 

ASSIGN Muff TO * 

MAT Avl= Data 
SUBEND 

SUB Filter(A(*),Fil) 

DIM E(1 : 1024) 

PRINT CHR$(12) 

PRINT USING ■////////’ 

PRINT ’ FILTERING THE SCATTERING PROFILE’ 

FOR J=1 TO Fil 

FOR 1=1 TO 1014 



B=I+1 


ITHIS IS A Synetric Hovinq Averaqe 


C=I+2 


! TYPE OF DIGITAL FILTER, EACH 


D=I+3 


'DATA POINT IS EQUALLY WEIGHTED 


F=I+4 


IIN THIS CASE BUT THIS CAN BE 


G=I+5 


ICHANGED IF ONE DETERMINES THAT 


H=I+6 


(FEWER POINTS WITH UNEQUAL WEIGHTS 


K=I+7 


IWOULD BE FASTER OR GIVE BETTER 


L=I+8 


'RESULTS, THIS TYPE OF FILTER WAS 


N=I+9 


IUSED SINCE IT INTRODUCES NO PHASE 


P=I+10 


ILAG (Angular Resolution) . 



E(G)=(A(I)+A(B)+A(C)+A(D)+A(F)+A(G)+A(H)+A(K)+A(L)+A(N)+A(P) )/l 1 
NEXT I 
MAT A= E 

NEXT J 
SUBEND 
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